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We consider ultracold quantum gases of scalar bosons, residing in a coupling strength-density 
regime in which they constitute a twofold fragmented condensate trapped in a single well. It is 
shown that the corresponding quantum states are, in the appropriate Fock space basis, identical to 
the photon cat states familiar in quantum optics, which correspond to superpositions of coherent 
states of the light field with a phase difference of tt. In marked distinction to photon cat states, the 
very existence of matter wave cat states however crucially depends on the many-body correlations of 
the constituent particles. We consequently establish that the quadratures of the effective “photons,” 
expressing the highly nonclassical nature of the macroscopic matter wave superposition state, can be 
experimentally accessed by measuring the density-density correlations of the interacting quantum 
gas. 


The eponymous cat quantum states have originally 
been constructed by Erwin Schrodinger to stimulate 
deeper thought on the sometimes surreal aspects of 
macroscopic quantum mechanics [ 1 ]. Present-day tech¬ 
nology allows for going beyond the pure gedankenexper- 
iment of Schrodinger’s day to create, within the realm of 
quantum optics, small cat states, Schrodinger “kittens” 
[2]. They consist typically either in the superposition of 
coherent states of light [3] with a phase difference of tt, 
coined photon cat states, or in states \N0) + \0N) of the 
NOON type [4]. Highly entangled photonic cat states 
constitute a possible basic building block in quantum in¬ 
formation architectures [IISHS]. 

We propose in what follows a novel species of 
Schrodinger cat states, which are in distinction to pho¬ 
tonic cat states relying on quantum many-body correla¬ 
tions, and hence owe their very existence to strong inter¬ 
actions of the constituent bosons. They are generated in 
a scalar (one component) gas of massive bosons trapped 
inside a single (harmonic or box) trap. The many-body 
states required to produce the cat states correspond to 
twofold fragmented condensates, for which the single¬ 
particle density matrix has two 0{N) eigenvalues [TO]. 
Condensate fragmentation m is by now firmly estab¬ 
lished as a many-body effect not describable within a 
mean-field theory of the classic Bogoliubov type [HUS]. 
Specifically, for sufficiently large and positive contact in¬ 
teraction coupling, it derives from the broken transla¬ 
tional symmetry and localization in a single trap m- 
Therefore, strongly interacting matter wave cat states 
(termed SIMCAS in what follows) can be experimen¬ 
tally accessed, starting from a Bose-Einstein condensate 
(BEC), by tuning the coupling with a Eeshbach resonance 

m- 

Below, we shall demonstrate that SIMCAS, in the ap¬ 
propriate Eock space basis, form many-body states es¬ 
sentially indistinguishable from photonic cat states, for 
sufficiently large (but still mesoscopic) number N of gas 
particles. The degree of fragmentation of condensate co¬ 
herence will be shown to be directly related to the de¬ 


gree of macroscopicity of the coherent superposition of 
many-body states, quantified by the superposition size 
of the quantum state pTUTO] , By evaluating quadra¬ 
tures [5], we show that the macroscopicity of the mat¬ 
ter wave quantum superposition is directly measurable 
through density-density correlations after time-of-fiight 
expansion (TOE). 

We emphasize that SIMCAS live in a single trap and 
not in a double well, and are obtained for a scalar, not 
a two-component gas, with repulsive interactions, which 
distinguishes them from previously suggested cat state 
implementations in dilute ultracold quantum gases, cf. 
[2M24] . Eurthermore, SIMCAS are the ground states of 
the trapped Bose gas for relatively moderate interaction 
couplings, at which the single BEC crosses over to a two¬ 
fold fragmented condensate [25] (the gas remaining suf¬ 
ficiently dilute to ensure that three-body recombination 
rates remain small). This renders SIMCAS different from 
the collapse-and-revival type quantum optics proposals 
to generate cat states in a Kerr medium cf., e.g., [26l|27|, 
or involving scattering of matter waves at barriers [28] . 

As a first step, we expand two-mode fragmented states 
in the Eock space basis of elementary bosons, 

N N 

\^) = j2ci\N-i,i), ^ia|2 = i, (1) 

/=0 1=0 


where the Eock space basis vectors \N — l,l) = 
two-mode Hamiltonian, with interaction part i7int = 


I vac), with I vac) the particle vacuum. A 


^alalaodo + ^a\a\aiai + ^ ^doagaiai + h.c.j + 

^a\aialao (assuming that A 3 G 3?), valid, e.g., for 
contact and dipolar two-body interactions when the two 
modes have even and odd parity, respectively, generically 
leads to a distribution of the Ci modulus which is, in the 
continuum limit, a Gaussian (of relative width oc I/\/]V), 
yielding fragmented condensate states for A 3 > 0 and 
Ai + A 2 + 2 A 3 - A 4 > 0 (see detailed discussion in M)- 
The distribution eenter Iq := {a\ai) G [0,N] quantifies 
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the degree of fragmentation = 1 — |Ao — Ai|/A^, where 
Xi are eigenvalues of the single-particle density matrix. 
We note that the exact functional form of the Ci distri¬ 
bution is not important for the following argument. The 
conditions are (a) 0 Iq N, (b) negligible weight of 
Cl at the boundaries I = 0, A/", and (c) small relative dis¬ 
tribution width. The fragmentation degree is = 2lo/N 
when lo < N/2 (assumed below without loss of general¬ 
ity) and T = 2(1 — lojN) when N/2 < Iq < N [141 [31]. 

Self-consistent solutions of the quantum many-body 
problem have indeed also found such two-mode frag¬ 
mented states for specific Hamiltonians, e.g. for quasi- 
one-dimensional (quasi-ID) BECs at sufficiently large 
contact interaction couplings [25l [29] . The large overlap 
of the modes in a single trap, crucially, leads to maroscop- 
ically large interaction-induced pair-exchange processes, 
— (dgagaidi+h.c.) ^ 0 {N‘^) which stabilize (rather than 
destabilize) SIMCAS, as shown in [30]. The existence 
of negative macroscopic pair coherence due to a pair- 
exchange coupling As ^ 0 {Ai, A 2 , A 4 /) is thus a distin¬ 
guishing feature of SIMCAS. 

To proceed, we construct the two ladder operators 

b = lim — , ^ dndi, b' = lim — , ^ ntdo, (2) 

where Ni = dta^. The ladder operators and their 
Hermitian conjugates b^,b'^ convert particles between 
the two macroscopically occupied modes, according 
to b\N-1,1) = Xl\N + , b^N-l,l) = 

\/l + 1\N — I—, I + 1), and similarly for b',b'^. The e 
regularization is introduced for the (finite N) singular¬ 
ity created when ¥ (6'’*') acts on |0,A^) (|A^, 0)), which is 
singular because there is no particle to be transferred to 
mode 1 (mode 0) for this state. 

We demonstrate in the following that a twofold frag¬ 
mented state can be written as the superposition of trun¬ 
cated coherent states of the ladder operators. These co¬ 
herent states live in a finite-dimensional Hilbert space 
[32] [ 33 ] for particle-number-conserved states of the form 
0. The ladder operators approximately satisfy bosonic 
commutation relations, due to 1 — (T| [6, 6’*'] |T) = (A^ + 
1 )\Cn\^ and 1 - (T| [6',6't] |T) = (Ar + l)|Cop, provided 
ICoUCivI «l/v^. 

The truncated coherent states \/3) of b are defined by 

b\(3) = b ^ 10 j = P \N) , (3) 

where |/) := \N — l^l). The normalization factor is given 

by A/s = exp(-|/3|V2)^ r(jv+i!|i3W ’ 
complete gamma function r(s,x) = e~^ dt. The 

state I/d) equals the usual coherent state in a infinite¬ 
dimensional Hilbert space when for A/" ^ oc the second 
term on the right-hand side of Q becomes negligible; in 


fact I/d) becomes very close to a proper coherent state 
already for moderate values of N [34] . Similar considera¬ 
tions hold for the truncated coherent states of the pair 6', 
I/d'). Hence both |/d) and |/d') represent a finite- A^ gener¬ 
alization of the concept of coherent states. The quasipar¬ 
ticles created by b^ or, as an equivalent choice (also see 
below), by 6'’*', will thus assume the role of the “photons.” 

The fragmented condensate states are transparently 
described in terms of phase space states [36] [37] , where 
the relevant phase is conjugate to the occupation number 
difference of the modes. For large N, fragmented conden¬ 
sates correspond to two phase space states oc 

N — /odg + I vac), separated in phase (j) by 

exactly tt m- We now show that individual phase-space 
basis elements are well approximated by truncated co¬ 
herent states of b (6'); we will then conclude that the 
fragmented ground state can be expressed as a superpo¬ 
sition of antipodal coherent states, Eq. 0 below. 

According to m, a two-mode many-body state 0 
can be expressed in terms of a phase-space basis as 

1^) = y* ^ 10, N, lo) , (4) 


where the relation of the phase space basis states cen¬ 
tered at lo to Fock space basis states is given by |/) = 
fo^ ^ 10, A^, /o): with the normalization fac- 

tor AfNM = 0Inverting the re¬ 
lation we have \(l),N,lo) = 
with Q, defining /d = |/d| exp[i0^]. 


|0,A^,/o) 



W 


10 ). ( 5 ) 


Now, by setting |00 = /q, that is by choos¬ 

ing the mean “photon” number to be equal to 
the (7/-distribution center, we obtain |0, A^,/q) = 

/o’ (E, I« ■ T1.0 

factor {N — lo)^~^/{N — /)! is proportional to a Pois¬ 
son distribution centered around N — lo^ well approx¬ 
imated by a normal distribution with both mean and 

variance N — lo. Then Yl^o \/ ^ w-z )\~ e~oc 

Ylf=o e~ This resembles a discrete 

Fourier transform of the normal distribution except that 
I > 0 and the summation range is finite. 

We thus find, with the choice |00 = /q? that a phases- 
space basis state maps to a truncated coherent state of 
the ladder operators, i.e., |0, A^,/q) — 10), as long as |0| 
and hence the fragmentation do not become too small, 
in which case the details of the Ci distribution at the 
boundary I = 0 matter and the concept of a 0-coherent 
state breaks down. An analogous argument can be ap¬ 
plied to 10') 10, A/", /o) [34]. Below, we restrict ourselves 
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to the truncated coherent states associated to the ladder 
operator 6, with analogous conclusions holding for b'. 

Using 1/3) /o) and |-^) + tt, /q), we 

thus make the ansatz C /3 |/d) + C-^ \ ~/3) for the 

fragmented condensate many-body state. To identify the 
coefficients, we observe that according to [HUSH], a frag¬ 
mented state can equivalently be written as a superposi¬ 
tion of two (for large N) degenerate many-body states 

|T) = c (|even) + uex.p[iOK] |odd)), (6) 


where |even) only contains even I coefficients Ci from 
the state 0 and I odd) only odd /, and |cp(l H- 
\u\‘^) = 1. By ^ and the fact that <pi = aig{Ci) = 
arg (C/ 3 -h C_/ 3 (-l)^) + Y (for 0/3 = 7r/2), it is es¬ 
tablished by straightforward algebra that we can write 

+ fur- 


even) 


\/ 2 (l+e- 2 |/ 3 | 2 ) 


eL % (i‘(i - (-1)')) 10- 


The overlap {—/3\/3) = ^ for 

large and we obtain even and odd superpositions 
of 10) and 1-0), |even) - (10) + |-0)), 

In summary, we have established that a fragmented 
condensate state can be written as a superposition of 
“photonic” truncated coherent states of the ladder oper¬ 
ators defined in 0, 




— ^0 5 
- 1/2 


(7) 


The 


where JV = (l-h -t-2r cos^exp[—2|00]) 

phase 0 and coefficient r in terms of u, Ok in 0 are found 

from rexp[i0] = ^ ^0 = \/111 -1% - 

In the limit |00 1, 


r = 


(1 — u^) + 2m cos 


(1 -h u‘^) + 2ixsin6>K 


0 = tan 


f 2 ucos0k 
1 — u‘^ 


.( 8 ) 


We thus find, in the equal-weight case r = 1, that the 
fragmented two-mode many-body states represent the 
analogue of photon cat states [SI]. For the latter, the 
truncated coherent state |0) of the massive bosons is re¬ 
placed by the coherent state of, e.g., a cavity photon field. 

It was shown in m that a stable variety of frag¬ 
mented two-mode many-body states is obtained for r = 

0 = 7rl2,37rl2. In the quantum optics literature on 
photon cat states, due to their potential for creating en¬ 
tangled coherent states, cf. the review [39], the states 
I even), |odd) are frequently studied, which correspond 
tor = 1, 0 = 0,7r. They are susceptible to perturba¬ 
tions when 100 1, which has been borne out for the 

presently studied ’’photonic” case in an interacting gas of 
massive bosons [30|. Finally, r 1 and r ^ 1 yield |0) 
and |—0), respectively, which represent nonfragmented 
condensation. 


The variances of effective position and momentum vari¬ 
ables (quadratures), in ladder space, are obtained from 


± A(S±Sf) 


4|/3|2cos^^^- hi, “-h 
4|/3|2sin^(/)^ (rf^) +f’ ““ 


for 100 1, and where A(b±b^); = {(b±b^)^) — (b±b^)^. 

The superposition 0 is a proper cat state for r ^ 1, 
with large quantum fluctuations of “position” and “mo¬ 
mentum.” In (|^, 0/3 can be considered as a parameter 
determining in which direction of “photonic” quadrature 
fluctuations are squeezed. 

The size of the SIMCAS (that is their degree of quan¬ 
tum mechanical superposition macroscopicity) is, in ac¬ 
cordance with the quadratures 0, determined by |00. 
The size of coherent state superpositions such as 0 is 
maximal when the overlap (—0|0) exp[—2|00] is min¬ 
imal [niiin], the size Ai being a simple polynomial func¬ 
tion of the overlap, AT ^ (1 — exp[—2|00])^. Thus, due 
to 100 =7’Ar/2, 

M ^ {1 — exp[—(r = 1). (10) 


The exponential dependence of the superposition size AT 
on T highlights the distinctly nonclassical character of 
a fragmented condensate in comparison to a BEG (sin¬ 
gle condensate). This will be manifest already for rel¬ 
atively small values of when N is mesoscopic, as in 
experimentally realized quantum gases. We remark that 
there exist alternative measures for assessing the macro¬ 
scopicity of superpositions, which rely on calculating the 
quantum Fisher information gnun], but obviously also 
involve the overlap (—0|0) in an essential manner. 

In marked distinction to superpositions of photon co¬ 
herent states, the superposition size AT strongly depends 
not only on a particle number, N^ but also, via on 
the strength of the microscopic (two-body) interactions 
in the constituent system dilute quantum gas; AT in fact 
vanishes exponentially fast for weakly interacting systems 
where IF ^ f). The macroscopic superposition size of the 
SIMCAS is therefore a genuine many-body effect |42|. 
We note in this regard that the macroscopicity AT for 
our matter wave cat is assessed by the rather straight¬ 
forward measurement of density-density correlations (see 
below), and not the elaborate quantum state tomography 
required for proper photonic cat states. 

Next we establish the relation of the “photonic” 
quadratures 0 to the density-density correlations in 
the interacting quantum gas. The nonclassical char¬ 
acter of the superposition 0, in conjunction with its 
many-body origin, will, by this means, become experi¬ 
mentally verifiable. We consider the density-density cor¬ 
relations in an effectively one-dimensional (ID) system 
(also see the below discussion on the experimental im¬ 
plementation), Ap 2 {z,%) := {p{z)p{z')) - pi{z)pi{z'), 
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where the density pi{z) = (plj^{z)%jj{z)) . We take the ap¬ 
proximation that (^p{z)p{z')^ 'Zi {pjj\z)%l)^{ z')%1){z')%Ij{z)')^ 
which holds true for (^i) ^ 1. For simplicity, we as¬ 
sume that the two macroscopically occupied field op¬ 
erator modes have even and odd parity, respectively, 
^{z) = '0o(^)^o +'0i(^)di, with '0o(^) = '0o(-^) and 
'0i(^) = and are real, ^i{z) G M. Bearing in 

mind that a fully self-consistent solution of the many- 
body Schrodinger equation will yield the true orbitals 
cf., e.g. [m |T3l Uni [43], for illustration purposes we 
take them to be the ground and first excited states of 
the harmonic oscillator; 'ipo{z) = 7r“^/^exp [—z‘^/2] and 
'ipi{z) = 7r~^^^\/2zexp [—^^/2], where the 2 : coordinate 
is assumed to be scaled by a suitable measure of (half 
the) extension of the cloud. 

For \P\‘^ ^ l,r = 1, we obtain the correlation func- 
tion Ap2{z,z') = i;o{z)'ipi{z')i^o{z')'ipi{z)[2{dla\aiao) + 
(aQfflQaiai) + (appoao) — ((fflofli +aiao))^]- Converting 
this to “photonic” ladder space, we get 


Ap 2 {z,z') ~ iJo{z)i’l{z')'lpo{z')'lpl{z) 


(6 + St))b - ((7 Vo'/2 (6 + b^))y 


( 11 ) 


The density-density correlations in the strongly interact¬ 
ing gas therefore map out the quadrature in the first 
line of ^ when the fluctuation of occupation numbers is 
small, i.e. for large N at finite T. 

Before TOF, 0/3 = |^, with no distinct correlation sig¬ 
nal present (the negative pair coherence of a fragmented 
state enforces sgn((7/+2C'/) = —1 (311 [38]. which in 

turn implies the condition 20/3 = tt + 2n7r, where n is an 
integer IS]). After expansion of the cloud, it has been 
shown in m that TOF of the cloud in the axial direc¬ 
tion effectively performs a rotation 0/3 ^ 0/3 — 7r/2 of the 
relative phase of the two modes, cf. the definition of b in 
(2) involving djai. This results in 0/3 = 0, and we con¬ 
clude from (9) and (11) that a characteristic correlation 
signal develops [44]. We plot the density-density corre¬ 
lations according to the state 0 in Fig.[^ which maps 
out with increasing accuracy the exact correlations of the 
underlying fragmented many-body state 0 when 10 p in¬ 
creases [ 34 ]. One clearly recognizes the emerging strong 
correlation signature of an increasingly macroscopic SIM- 
CAS for larger “photon” number |00: The larger the de¬ 
gree of fragmentation, the bigger the cat becomes. To 
identify the ‘dead’ and ‘alive’ parts of the superposition 
state, we plot in the bottom part of Fig.the density- 
density correlator (^p{z)p{z')')^ for three limiting cases of 
superposition weights in 0, r = 0, 1 , 0 c, yielding |0), 
^(10) +exp[i^] |—0)), and |—0), respectively. 

For the experimental realization of SIMCAS, quasi- 
ID systems [45] are favorable, since they possess larger 
degrees of fragmentation when compared to higher¬ 
dimensional gases m- SIMCAS exist in a interaction 



{p{z)p{z')) 


ol^|o 


-1 0 1 z -1 0 1 z -1 0 1 z ° 

\-f3) (‘Dead’) |/3» |/3} (‘Alive’) 

FIG. 1. Emergence of a “photonic” cat state superposi¬ 
tion in an ultracold quantum gas containing A = 100 mas¬ 
sive, interacting bosons. Top: ID density-density correlations 
Ap 2 {z,z'), scaled by where Z is the extension of the 

cloud after axial TOF (Z oc t in the long time limit [31]), 
which is also the scaling unit of z,z'. From left to right, the 
fragmentation degree increases according to = 0.1, 0.2, 0.4 
(100 = 5,10, 20, with r = 1). Bottom: Density-density corre¬ 
lator (p{z)p{z')) (in the same N^jZ^ units) for |00 = 20, and 
r ^ 0 (left), r — 1 (middle), r ^ 00 (right). The correlator 
does not depend on 0 in the large |00 limit. 


coupling-density region of the phase diagram relatively 
close to the BEC domain and far from the extreme Tonks- 
Girardeau regime of ultralow density and very large cou¬ 
pling gaiiT] (in which the number of fragments equals 
the number of particles [25]). The requirements on gas 
density and contact interaction coupling constant g to re¬ 
alize SIMCAS are therefore more moderate than for the 
Tonks-Girardeau gas, and should be readily accessible in 
experiment. Finally, the phase Ok and amplitude r in the 
superposition 0 can be engineered by a rapid sweep of 
interaction couplings, as demonstrated for Ok in [38] . 

Our analysis reveals the highly nonclassical character 
of fragmented condensate many-body states, as opposed 
to the essentially classical BEC contained in the same 
trap. We therefore anticipate applications of SIMCAS, 
inter alia^ in quantum metrology [48] , where the Cramer- 
Rao bound furnishes rigorous bounds on the precision to 
which parameters in the Hamiltonian can be measured 
[49] . The quantum superposition macroscopicity of SIM¬ 
CAS manifest in the density-density correlations of the 
bosonic gas thus establishes a potential many-body re¬ 
source of parameter estimation theory. 

We thank Ty Volkoff and Daniel Braun for helpful dis¬ 
cussions. This research was supported by the NRF Ko¬ 
rea, Grant No. 2014R1A2A2A01006535. 
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SUPPLEMENTAL MATERIAL 

We here provide a more detailed analysis of both the accuracy of using \p) as approximate (truncated) coherent 
states, as well as the accuracy of the identification of phase-space basis states with them, |0, A, /q) — |/^)- Furthermore, 
a more extended description of the TOF evolution of density-density correlations is presented, and the accuracy of 
the cat state ansatz |4^) = M{\P) + |~/^)) for fho fragmented condensate many-body state is assessed. 


Accuracy of Treating |/ 3 ) as a Coherent State 


In the main text, we introduced in Eq. (2) ladder operators, which represent approximate bosonic annihilation 
operators, as follows 


t) = lim — —— 

VNo + e 


al&i, 


b' = lim — ^ 

V A + e 




and the corresponding truncated coherent states as 


A /3' 


JV 




/?' 


■fN-l 


1=0 


\N-l,l), 


(SI) 


(S2) 


toVW^- 

where = exp(-|/3|V2)^ r(jvA7^|'-^r Truncated coherent states were previously treated, e.g., by [32]. They are 
finite-dimensional-Hilbert-space versions of the bosonic annihilation operators d and coherent state ja) in quantum 
optics. Here, we aim at finding a quantity which assesses the accuracy of treating b and \P) {b' and |/3')) as bosonic 
annihilation operators and their corresponding eigenstates, 

[6,St] = l, b\p) = p\p), (S3) 

and which allows us to numerically evaluate the precision to which the above two relations hold. 

We assume, without loss of generality, (djdi) < A/2 with |/3p < A/2 because, as will be shown in the following, 
more robust sets of bosonic operator and their eigenstate are {6, \^)} when (djdi) < A/2, and {6', \P')} when 
(djdi) > A/2. In addition, we note that there is always the freedom to choose 6, \f3) or 6', \/3') to describe a given two 
mode system. With respect to an arbitrary (normalized) two-mode state |4^) 




1=0 


1=0 




(S4) 


1=0 
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I /31 

FIG. SI. Plots of —(indicated by loss of \Ci\^ on the vertical axis), with |/3|^ = iV/4 (blue) and |/3|^ = N/2 
(black), for N = 2b (left) and iV = 50 (right). 


the accuracy of (S3) can be assessed by evaluating the following two quantities 


l-(^|[fe,St]|$) = (7V + l)|CivP, 


((/?|/3*-(/?|St) (/3|/?)-b|/3)) 


\N 




= Ai 


m 


(S5) 


It is readily observed that the accuracy of 6, as a proper set of bosonic operators depends on the state |^) which is 
considered. In the main text a twofold fragmented state is discussed with the ansatz |^) = N'{\I3) \ (also see 

below), and since the \Cn\‘^ distributions of \/3) and \—/3) are identical, it is sufficient to consider |^) = \/3). Eq. (S5) 
becomes 


1-{P\0-b^b)\p) = {N + l)Ay 


| 2 iv ((^1- (^1 6t) (/3 1^) - 6 


m 


l/?P (/?!/?) 


_ ,2 \pr 


(S6) 


Higher order expressions in such as bWb^ or the repetitive action of b against |/3), can be similarly assessed. 
They induce the “loss” of certain amplitudes Cw-i, * * • , Cw-n+i from the state vector where 


N 

r I/?) - {br 1/3) = 13^Ap Y. 

l=N—n-\-l 


Vv. 


\N-l,l). 


(S7) 


Therefore the robustness of the coherent state representation can be defined by the degree to which a state of interest 
|\l/) is not deformed from the loss of the amplitudes Cjv, • • • , C^-n+i due to a n-times repeated action of b on it. 
Specifically, for |\l/) = |/3) this robustness can be quantified through 


AT 

Z-^l=N—n-\-l 


\lf _T{N + 1, |/3|2) - T{N + 


A/" 

Z^/=0 


/! 


r(iv + i,|/3|2) 


- exp(-|/3|" 


N 

' ^ n 

l=N—n-\-l 


21 


(S8) 


which represents the lost fraction of the Ci due to repetitive action of b. We used here r(A^ + l,|;5p) r(A^ + l) 
which holds within 1% error for N > 10 with given \/3\‘^ < N/2. In Fig. |Sl[ we show plots of the quantity defined 
in (S8) for different A/", n, 1/31“^/N. We see that as N gets larger and 1/31“^/N gets smaller, the robustness significantly 
increases. Specifically, for A/" = 25, = O.bN there is already a negligible loss of Ci amplitudes from the state even 

upon n = 5 times acting with 5 on it. 


Comparison of |0,iV,/o),|^),|/3') 


We argue in the main text that 10, A^,/q) — |0) for small Iq = |/3p, and |0, A^,/q) — \P') for small N — Iq = \f3'\‘^ 
with 0 = 0/3 and 0 = —0/3', cf. the discussion after Eq. (5) in the main text and Eq. (S9) below. We present in Eig.|S2| 































FIG. S2. Plots of the \Ci\^ distribution of |0, A/", /o) (red), |/3) (blue), \I3') (green) with N — 100 for different Iq — |/3|^ = 
N — l/3'l^ = O.liV, 0.2iV, 0.3iV (Top from left) and h = = iV — |/3'|^ = 0.7iV, O.SiV, 0.9A^ (Bottom from left). 


the distribution for 10, A^,/q) , |0), \P')^ respectively, with different /q = in the case N = 100. We conclude 
that 10, A^,/o) — 10) is confirmed for small /q = 100, and |0, A^,/q) — |0O small N — Iq = |0'p with 0 = 0^3 and 
0 = -0/3'. 

Furthermore, we have verified that in the large N limit, the truncated coherent states to increasingly good accuracy 
represent the corresponding phase-space basis vectors even when lo ^ limit of maximal fragmentation, 

T ^ In more detail, as |00 approaches |00 = A^/2 [T = 1), the |C/0 distribution of the phase state |0, A^,/q) 
becomes wider than that of |0) or |0'). However, |0, A^,/q) — |0) still holds. From equation (5) of the main text, 
which transforms |0, A^,/q) into the |0) basis, we have 


|0,A^0o) = 


m 


A 


^ d<Pi3 ( 
i/3 Jo V 7 


{N-l)\ 


p27T 

Jo 


27T 


Cs 


(S9) 


S3 


with |/3|2 = la where = 0^^ (E; 

distribution, is plotted for Iq = |0p = A/'/2, and with various values of 

with increasing N, even when A^ = 25, a rapid convergence of | to one peak at 0/3 = 0 occurs. 


the modulus of the coherent state phase 
particle number N. It is apparent that. 


Time-Of-Flight Expansion (TOF) and Phase Rotation 

We present here a more detailed analysis of the phase rotation incurred by TOF. 

After TOF, the assumed even-odd parity of the modes is preserved, 'iPq{zA) = = “'0i(—^ 5 ^)* 

Defining the phase variable cp such that 

:='tpi{z,t), ipQ{z,t)'ipi{z,t) eM. (SIO) 

the wave functions ‘ipo{z,t) and ‘ipi{z,t) share a common phase factor [35], keeping (ata^y (atatafca;) {i,j, k,l = 0, 1) 
invariant under time evolution in the noninteracting limit of TOF. Before TOF (t = 0) (p = 0 and after TOF (large 
t) (p = —1^12 [3T] . 

Note that the value of p could be both p ox p P n since times a real number is again real number. We write 
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FIG. S3. Plots of the distribution of 10, A/",/o) (cf. Eq. (5) in the main text and the discussion that follows it), for 

/o = 101^ = N/2 with N = 25, 50,100. 


down the density-density correlation function Ap 2 {z^ z'^t) as 

Ap 2 {z,z',t) = \ipo{z,t)f\ipoiz',t)f 0ajajaoao) - (aJao)(aJao)) + \ipi{z,t)f\ipi{z',t)f 0a|a|aiai) - (a|ai)(a|ai)) 

+ i\i’o{z,t)f\i^iiz',t)f + |V’i(2,t)nV’o(^',i)H ((ajabi^o) - (aoao)(4«i)) 

+ {\'ipo{z,t)\^'ipUz'+ \'<Po{z',t)\^'ipo{z,t)'ipi{z,t)) ((ajajaoai) - (aJ)ao)(4«i) 


+ {\'ipi{z,t)fi>Q{z',t)'ipi{z',t) + \i>i{z',t)fi>Q{z,t)'ipi{z,t)) e*‘^ 0aja|aiai) - (ajai)(a|ai)) + h.c. 

+ ^pQ{z,t)'tpi{z,t)^pQ{z',t)^pi{z',t) 2(ala[aiao) + (alalaiCLi) + e“^*‘^(a|a|aoao) - + e“*‘^a|ao)) 


(Sll) 


Furthermore, we can take '0i(z,t) as the new '0i(z,t), yielding 2 po{z,t) 2 pi{z,t) G M, by changing the many-body state 
1^) as follows 


N N 

|V|/) =J2Ci\N-l,l)^ = J2 \N -1,1). 

z=o /=0 


(S12) 


This rotation of the state in (S12) corresponds to a rotation of the phase (pp 


Ui UiiC CippiUT^iliidUC 


1/3) 


(S13) 


For a general twofold fragmented state we then get 
1 


vT 


( 1 / 3 ) 


A9 


-/3)) 


VT 


(|/3e*‘^) +re*^ \-/3TV) ■ 


(S14) 


Since (pp = it/ 2 before TOF, it is rotated to 0/3 = 0 after TOF. We can therefore summarize as follows: TOF evolution 
is equivalent to (p = 37r/2 = —7r/2 rotation of state defined in ( |S12[ ), and we can calculate Ap 2 {z^z') from (Sll) for 
both before TOF {p = 0 or (j)p = 7t/2) and after TOF {p = —tt/2 or 0/3 = 0). 


Validity of the |T) = ^(1/3) -h |—0)) Ansatz 


As a main ingredient of our analysis, we propose an ansatz, |T) = A/'(|/3)+r*^ 1“/^))? describe a twofold fragmented 
state (Eq. (7) in the main text). To test the validity of the ansatz in more detail, let us first specify the Ci distribution 
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of a general twofold fragmented state. In m, it was shown that \Ci\ has a Gaussian distribution of mean Iq and 
variance with a ^ 0{-s/N) as long as the continuum limit for the Ci distribution holds. Furthermore, a general 
twofold fragmented state can also be written as [38] 

1^) = c (|even) + i^exp['i6>K] |odd)), |c|^(l + = 1, (S15) 


which leads to 


\Ci\ = V2|c| (2^/2)v^ ^ ' ^ 

IQ I = \/2|c||M| (2^^2^i/4 e~* 4 ,.°^ e*'^' for odd I 


(S16) 


where Ci = \Ci\e^^^ ^ and (pi obeys (pi -\-2 = + tt, (pi -\-4 = pi from the condition sgn(C/(7/+2) = —1 [Hj- We therefore 

have 


pi^i = pi^Ok for even I 
pi^i = pi ^ TT - Ok for odd I 


(S17) 


We now examine to which extent using |^) = PPdP) + |~/^)) is accurate for the calculation of the density-density 


correlations Ap 2 {z, z') by calculating (Sll) with ( |S16[ ) in terms of a, Iq and c^u. We then compare the corresponding 
result for the density-density correlations obtained by using the exact large N fragmented state with what we derive 
using Eq. (11) and the quadrature fluctuations of Eq. (9) in the main text. As we will show, we obtain two results 
whose difference depends on the value of the Gaussian width a. 

To calculate ( |S11[ ), we perform the following approximations. When the Gaussian \Ci\ distribution is well localized 
in the interval [0, N] and N is large enough to approximate Ci to be continuous, then it is permissible to approximate 
the expectation value of f{N — /, /), which is a polynomial of y/N — P \//, as 


I z=o,2,--- yzira 




"f 


f{N-l,l) 


1 


\/27rcr^ 


TTcr^ 

20-2 


+ 2|c|>|2 ^ nN-l,l) 




V 2 


TTCr^ 


(S18) 


dl. 


The quantities (dQdo),(d|di),(aQdQdodo),(a|d|aidi) and (dQdjdido) are directly determined from the above expres¬ 
sion. We have 


(djdjdodo) = {N - Iq)^ - {N - Iq) + cr^, (a|d|didi) = - /q + cr^, (dJd|diao) = {N - lo)lo - 

with (agdo) = N — /o,(d|di) = Iq. 

The quantities (aQdi),(dQdQaodi) and (dQdjaidi), in turn, are related to the sum 


a 


(S19) 


^ /(TV - /, /)cra+i ^ 2|cp|^| ^ /(TV - /, (S20) 

i i 

the last relation holding when \Ci \ IC'z+i I • Under the proviso that pi -\-4 = pi and \Ci \ slowly vary so that \Ci\ |Cz+41^ 

we can carry out the summation over edT>z+i-T>d separately from the I summation as 


^/(AT - ~ i ^/(fV - 

I ^1=0 I 

Therefore we have 


Ym-i,i)c:ci+, 


4 ^ 

z=o 








V 2 ' 


(i-ipr „ 
2ct2 dl. 


TTCr^ 


(S21) 


(S22) 


Erom the following string of phase values, pi = • • • , 0, 6>/e, tt, 6>/e + tt, 0, 6>/c, tt, + tt, • • •, where = 0 is chosen, we 
obtain 


1 

= zsin6»fe. (S23) 

^ 1=0 
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The function — now includes square roots of — / or /. It is therefore nontrivial to write down expectation 
values in terms of N^Iq^cf in closed form, with the exception of (djai) = 2\c\^uiy^ {N — lo)lo sinOk- Now, (dQaQaodi) 
and (ajdldiai) read 


/ oo 1 / 111 

{N -I - -)y' {N -1+ -){l + -) 


2 V27rcr^ 


(i-ior „ 
e 2 ct 2 dl, 


/ oo 2 I ^ 2 

(^ - 2 ) Y ^ + 2^71^ 

From the following integrals 

/ / , e 2 ct 2 dl = Iq. / 

7-00 \/w 7- 

we can power expand as 


:7rc7^ 


:e 2 <t2 dl. 


(S24) 


C» -I 

2 


. (lulo). 


\/2'kg^ 


e 2^2 (i/ = /g 7 - 


(djdjdodi) ~ 2|c|^msin6'fc (A^ - ?o)7(-^ “ ^ 0)^0 + 0(A'' - /q) + 0{a‘^) 
(dod|didi) ~ 2|c|^msin6'fe Iq^/{N - Iq)Iq + 0{N - Iq) + O(cr^) 


(S25) 


(S26) 


The first, highest order, terms are matching (aQao)(aQdi) oc {N — lo)y/{N — lo)lo and (aQai)(api) oc loViN-io)lo; 
we assume Iq < N/2. Finally, for (ao®o®i®i) ((oi«ho«o) is given by the complex conjugate), with slowly varying 
IC/ 1 , we have 


1 ^ pOO _ -I n_i \2 

(djajiaiai) - jYl / 7(Ar - - ^ + 1)^(^ + 1) ^—- e ^ dl. 

4 ^ J-oo v2iTa^ 

A expansion similar to that leading from ( |S24 ) to ( |S26| ), and \ Y^=o = —1, yields 

(djujaiui) ~ -(AT - lo)lo + 0{N - lo) + 0{a^). 

Summing up, Eq. ( |S11| ) becomes 

Ap 2 {z,z',t) = \lpo{z,t)f\lpo{z' ,1)1"^ (cr^ -{N - lo)) + \lpi{z,t)f\lpi{z',t)f (cr^ - Iq) 

+ {\‘4’o{z,t)\‘^\tpi{z',t)\‘^ + \xpi{z,t)\‘^\tpo{z',t)f) 

+ 2|cpMisin0fe {\'ipo{z,t)\^'ipoiz',t)'ipi{z',t) + {ipoiz',t)\^'ipo{z,t)'ipi{z,t)) [0{N - Iq) + 0{a^)] 

+ 2|c|^msin6»fe {\'ipi{z,t)\‘^ipo{z',t)'ipi{z',t) + \i)i{z',t)\‘^ipoiz,t)'ipi{z,t)) [0{N - Iq) + 0(cr^)] 

+ 4'ipo{z, t)'ii!i{z, 1)4)0{z',t)'ii!i{z',t) [(-^ “ ^o)^osin^ (l - (2|cpusin6»fe)^) + 0{N - Iq) + 0(a^)] 


(S27) 


(S28) 


(S29) 


Provided the orbitals 2 po{z,t) and have large overlap, as to be expected in a single trap, the following four 

terms have similar order of magnitude, 

4*o{z,t)4i{z,t)4*o{z',t)4i{z',t) - \4o{z,t)?\4o{z',t)\^, \4i{z,t)\^\4i{z',t)\^, \4o{z,t)?\4i{z',t)?- (S30) 

This leads to 

Ap 2 {z,z\t) = A4l{z,t)4i{z,t)4l{z',t)4i{z',t){N - lo)^osin^</? (1 - (2|c|^Msin0fe)^) + 0{N - Iq) + 0(a-^). (S31) 

Therefore, since the first term of ( |S31| ) is proportional to {N — lo)lo, and because we assume {N — lo)lo ^ A — Iq, {N — 
lo)lo ^ we obtain Ap 2 {z,z') to lowest order as 


Ap 2 {z,z') ~ 4:ipo{z,t)i’i{z,t)4o{z',t)4i{z',t){N - lo)losm^ ‘p{^- (2|c|^usin6»fe)^). 


(S32) 


To compare with the result in Eq. (11) of the main text, we employ the relation between r^O and specified for 
large |/d| by Eq. ( 8 ). Then 2 |cpi 4 sin^/c = 2u/{l v?) sin Ok can be expressed in terms of u^Ok as 


1 - 
1 + 


2uXp (2u 


1 + u^ 




(S33) 
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In the A /3 = y^(l + e“^l^l^)/(l — ^ 1 limit (equivalent to e ^ 0 and small overlap of |/3) and |—/d)), this 

becomes 


1 — 
1 + 


‘2iU 

— —5- sin^fe = 2|cpMsin0fe. 


(S34) 


The result does not depend on 0. Thus we see that when the overlap between |/3) and |—/3) gets smaller, the effect of 
0 on Ap 2 {z, z') becomes negligible. In this limit (S32) in terms of r, with \/3\‘^ = is given as 


Ap 2 (z, z') ~ {z, t)'ipi{z, t)ipo {z', t)ipi {z', t){N - lo)lo sin^ p 


4^2 


(1 + r2)2’ 


(S35) 


which agrees with the r = 1 result of the main text in Eq. (11) as long as {N — Io)Iq ^ N — Iq, . 

We conclude that {N — lo)lo ^ A^ —/q, and the large single-trap overlap between the orbitals '0o(^5 and '^ 1 ( 2 ;, t) 
are the conditions to utilize the representation |T) = A/’d/d) + \ —/3)) for a given twofold fragmented state. 








